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I. INTRODUCTION 

This is one of several articles describing theoretical 
calculations of first-forbidden [3 decay observables and 
other related weak-interaction variables with as high ac- 
curacy as is currently possible. Recent results have been 
reported for A = 50 §, and the A=133-134 @, and 
A=205-212 H regions. Here we report on the four A = 
11-16 decays of Fig. |l| and Table |. The main motive 
is to add understanding of the very large enhancement 
over the impulse approximation of the rank-zero (R0) 
axial-charge matrix element Mq observed especially in 
the lead region [M. 

The four (3 decays of Fig. |l| are the known decays in 
light (A < 37) nuclei that are fast enough to provide 
potentially reliable information on the medium enhance- 
ment of Mq . In addition to these four decays, an impor- 
tant part of this study is a consideration of the inverse of 

16 N(CT) JZ^ 16 O(0+), namely n~ capture on 16 lead- 
ing to the 0~ first-excited state of 16 N with — for the 
first time — a calculational precision comparable to that 
routinely used in /3-decay studies |llj . 

The light nuclei have always been and will continue 
to be the premier testing ground for our views on the 
structure of nuclei. Our main emphasis here will be on 
A = 16 nuclei. 16 has a fascinating and complex struc- 
ture since QTiuj, 2Tiuj, and Ahto excitations are manifestly 
apparent amongst the low-lying levels. The (0+2+4)?kj 
model of Brown and Green @ was an early, success- 
ful, and important description of these states. Recently, 
we described a large-basis shell-model diagonalization of 
16 in a (0+2+4) Tiw basis and 16 N in a (l+3)?iw basis 
Jl3|. Similar wave functions will be used in the 16 N(0~) 
«-> 16 O(0 + ) (3~ and [l~ calculations. The /3~ decays of 
16 C, 15 C and 11 Be will be treated in the more truncated 
(l+3)7kt> (0+2)^w space (A = 11,15) or the (2+4)Hlo 
— > (1+3)Huj space ( 16 C) with the effects of 4Uuj or bfiw 
added perturbatively. 

In Sec. II we give a review of the shell-model interac- 
tions used in the present study. In Sees. III. A and III.B 
we describe the calculation of the one-body (impulse ap- 
proximation) matrix elements, which enter in these rank- 
zero processes, and how they are combined to give the- 
oretical rates. The two-body (meson exchange) matrix 
elements are considered in Sec. III.C. In Sec. IV we 
describe the 16 N(0 _ ) <-* 16 O(0 + ) transitions and in Sec. 
V we take up the (3~ decays of 15 C, n Be, and 16 C, Fi- 
nally, in Sec. VI we discuss the conclusions concerning 
the enhancement factors for /3 _ and fi~ capture. 



II. THE WAVE FUNCTIONS 

Shell-model calculations were performed with the shell- 
model code OXBASH @. With OXBASH, spurious 
center-of-mass motion is removed by the usual method 
jl5| of adding a center-of-mass Hamiltonian H cm to the 



interaction. The shell-model studies start with the re- 
cently constructed interactions of Warburton and Brown 
[PH , which are based on interactions for the OplsOd shells 
determined by a least-squares fit to 216 energy levels in 
the A = 10-22 region assuming no mixing of nhuj and 
(n+2)hoj configurations. The OplsOd part of the inter- 
action cited in Ref. jl6| as WBP results from a fit to 
two-body matrix elements (TBME) and single-particle 
energies (SPE) for the p-shell and a potential represen- 
tation of the Op-(lsOd) cross-shell interaction. The WBP 
model space was expanded to include the 0s and Oflp ma- 
jor shells by adding the appropriate Oflp and cross-shell 
IsOd-Oflp two-body matrix elements of the WBMB inter- 
action Jrij ] and all the other necessary matrix elements 
from the bare G matrix potential of Hosaka, Kubo, and 
Toki 0. The Osjtf, and lp SPE were determined as 
described in Ref. [[L6| . Thus the WBP interaction is con- 
structed in a similar manner to the Millencr four-shell in- 
teraction described in Ref. jl9| , but reproduces the bind- 
ing energies of low-lying >l?ku levels in the A = 16 region 
with 2-3 times greater accuracy. 

Unless a complete model space is used for a given di- 
agonalization, the wave functions can contain spurious 
components. For a (0+2)?ia; calculation in 16 0, the first 
four oscillator shells comprise a complete basis, while 
for a (0+2+4)huj calculation the first six shells must be 
included for completeness (l^j2^]. We began our shell- 
model studies by diagonalizing the (0+2+4)?kj + T = 
and (l+3)7ia> _ -3 - T = 1 states of 16 in model 
spaces comprising both the first six and first four oscilla- 
tor shells. As discussed in Ref. Q, negligible difference 
was found in the wave functions and observables of in- 
terest between the calculations within these two model 
spaces. Thus the calculations reported here were per- 
formed in the four-shell model spaces — we emphasize 
that the results in six-shell model spaces would be essen- 
tially identical. 

In Ref. jl3| several methods for constructing the mixed 
Huj spectra were explored, each of which gave a reasonable 
level scheme for A=16. Since that work we have further 
explored the sensitivity of other observables such as the 
electron scattering form factors from the ground state of 
16 to + , 2 + , and 4 + excited states and the Ml decay of 
low-lying 1 + T = 1 states [^lj. We found that the method 
of lowering just the Atiu configurations for the positive- 
parity states tends to produce too much mixing of the 
4hui components into those states around 12-18 MeV, 
which are known — from experimental observables — to 
be mostly 2tku. The gap method turns out to be better in 
this regard because it simultaneously lowers the 2htu con- 
figurations with respect to the 4Tllu configurations. The 
needed reduction of the Op— IsOd gap — initially A(0p- 
lsOd) = 11632 keV — for the WBP interaction is 3.0 
MeV, which means that the 2hio states are shifted down 
by 6.0 MeV and the %hio states are shifted down by 12.0 
MeV. It is worth noting, as discussed in Ref. Jll|, that 
the amount that the 4hu> configurations need to be low- 
ered is essentially equal to the shift of the ground state in 
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going from the Ohu space to the (0+2+4)7iw space. Thus 
the view can be taken that the shifts are not free param- 
eters but with a given frcu truncation can be determined 
self-consistently. When the negative-parity T = 1 states 
are treated in a (1+3)7kj space it is perhaps reasonable to 
shift the 3huj configurations down by 12.0 MeV and the 
Yfou) configurations down by 6.0 MeV, and doing so gives 
excellent agreement with experiment for the energies of 
these states relative to the positive-parity T = states. 
Note that this is not the same as the pure gap method, 
which would involve shifts of 9.0 and 3.0 MeV for the 
2>hu) and ITwo configurations, respectively. Nevertheless 
we will refer to this as the gap method since the struc- 
ture of the negative-parity states are exactly the same 
in either case (since the structure only depends upon the 
6.0-MeV shift difference between the 3hu> and Vhu config- 
urations). We refer to the results obtained with the fitted 
potential for all matrix elements and the gap method as 
WBP. 

In addition, we have further explored the interaction 
y2hu fosfc mixes nhoj and (n+2)hu> configurations. In 
Rcf. [[D| this interaction (40 Op-lsOd TBME) was gener- 
ated from the potential obtained from the fit described 
above (WBP). However, since the core-polarization cor- 
rections may be different for Ifiio and 2tuv matrix ele- 
ments, the use of a common potential may not be en- 
tirely correct. As an alternative we have explored |2]J] 
the use of the Bonn G matrix 22 for the mixing inter- 
action V 2huJ . We find that the bare G matrix cannot re- 
produce the energy-level spectrum 16 with any amount 
of 2hu and/or Afiui shift, but that the interaction was ac- 
ceptable in this regard if it was renormalized by a factor 
of 0.8. Such a renormalization seems reasonable because 
the core-polarization corrections are found to reduce the 
average V 2huJ matrix elements of the bare G matrix by 
about this amount |23|] . We will refer to results obtained 
with V 2huJ replaced by 0.8 times the Bonn G matrix as 
WBN. For this interaction the needed reduction in the 
Op-lsOd gap was found to be 3.5 MeV so the the 2hu) 
and ATiuj components are lowered by 7 and 14 MeV, re- 
spectively. As for the WBP interaction, the spirit of the 
gap method was retained for the odd-parity states in that 
the difference between the STilo and IUlo shifts was set 
at 7 MeV. However the absolute shifts were arbitrarily 
set so as to minimize the difference in the experimental 
and theoretical excitation energies of the low-lying T = 
1 quartet. This, of course, has no influence on the wave 
functions of the states. 

The main reason for exploring the use of the Bonn 
y2huj interaction was that the WBP interaction did not 
provide a satisfactory description of the electron scatter- 
ing form factors and of the Ml decay of the low- lying 
1 + T = 1 states. The WBN interaction does consider- 
ably better in this regard and, in fact, provides the most 
satisfactory description of these phenomena of any re- 
alistic interaction that we are aware of Thus the 
results given here for the WBN interaction are the pre- 
ferred ones. The WBP results are given to provide a 



measure of the sensitivity to the interaction used. 



III. THE WEAK INTERACTION PROCESSES 

A. The one-body contribution 

Our concern here is with the one-body (impulse ap- 
proximation) matrix elements of the rank-zero (R0) axial 
current. These are the matrix elements of the R0 member 
of the spin-dipole operator and of the helicity operator 
75, which is commonly called the time-like component of 
the axial current, or the axial charge. With 75 replaced 
by its non-relativistic limit (good to order 1/Mn) the 
single-particle matrix elements (with the relative phase 
appropriate to (3 decay) are |2J-|2j 

Mg(jdf) = gAVz(jf\\MCi,<T] T\\\ji)O_ 



where 

Ct.i = 

and 
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and where Xce — the electron Compton wavelength di- 
vided by 2ir — is incorporated into Eq. (lb) so that both 
matrix elements have the dimensions of fm. The ji and 
j f are a short-hand notation for all quantum numbers 
needed to label the single-particle states. In Eq. (||), the 
Ji and T dependencies result from the reduced matrix 
elements in a spin and isospin space and \[2 is from the 
definition of the isospin operator. It is important to keep 
in mind that the matrix elements of Eq. (|l]) are her- 
mitian conjugates, which for harmonic-oscillator radial 
wave functions results in the identity 



E„ 



m e c* 



for HO. 
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The usual shell-model procedure is followed of combining 
these single-particle matrix elements M^(jijf) with one- 
body transition densities Dl (jijf) via 

m% = ]T M a R (jd f ) = J2 d r UdfWZUdf) (5) 



where the subscript R denotes the rank of the operator 
(R = in this study). The one-body transition densities 
given by 



(JfT f \\\[al®a Jf _ 



ATAJI 



\JiTi) 



[(2AJ+ 1)(2AT+ 1)]2 



(G) 
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contain all the information on the initial and final many- 
body wave functions. The J{Fi and JfTf are a short- 
hand notations for all quantum numbers needed to de- 
scribe the many-body wave functions and AT,AJ are 
multipolarities of the one-body operator, which in our 
case has AT = 1 and A J = R = 0. 

In the previous calculations for heavier nuclei J^-§] , the 
single-particle matrix elements were augmented by mul- 
tiplicative renormalization factors q a (jijf) that repre- 
sented first-order core-polarization effects. In the present 
calculation the model space is large enough to include 
all possible first-order core-polarizations, and so the one- 
body operators appropriate to bare nucleons are used. 

The Mft(jijf) and also the the two-body matrix ele- 
ments described in Sec. III.C are calculated with a com- 
bination of harmonic-oscillator (HO) and Woods-Saxon 
(WS) wave functions. The WS parameters were de- 
termined by extrapolation of results obtained by least- 
squares fits to nuclear charge distributions for 12 C and 
16 O p8| . In this procedure the separation energies were 
fixed as described below and the orbit occupancies were 
taken from (0+2+4)?ki; wave functions for 16 |L3| and 
12 C H]. For HO wave functions, we first found values 
of Tiuj for 12 C and 16 O that reproduced the root-mean- 
square charge radii, (r 2 )a , with these orbit occupancies. 
For 16 0, with (r 2 )^ = 2.730 fm @, the result is hto = 
13.60 MeV as opposed to 13.16 MeV, which is obtained 
assuming an 16 closed shell. Values of fiu) for A = 11 
and 15 were obtained from the A = 12 and 16 values 
assuming a linear dependence on A. 

The WS results depend on the separation energies S(n) 
for the (3~ parent state and S{p) for the daughter. These 
are related by 

S(p) - S(n) =Q/3 - 0.782 MeV. (7) 

where Qp — 0.782 MeV is the difference in binding ener- 
gies of the initial and final states. The separation ener- 
gies for a particular common parent state with excitation 
energy E x in the (A— 1,Z— 1) core are then 

S{n) = AE b (n) + E x - E t (8a) 



is(ls0d) — > 7r(0p); i.e., all others have large effective val- 
ues of E x and thus of S(n) and S(p). 

The method developed by Millener |30| and adopted 
to a similar calculation |l(| for 206 T1 — >^ 6 Pb was used 
to estimate the effective value of E x , (E x (j)), to use in 
Eq. (^). This method involves an inclusive calculation 
of the spectroscopic amplitudes for neutron pickup from 
the initial state and proton pickup from the final state 
to all possible core states and an evaluation of the effec- 
tive excitation energy of the core states from a consider- 
ation of these amplitudes and the resulting dependence 
of the matrix elements on E x . In this determination all 
core-state excitation energies greater than 10 MeV were 
fixed at 10 MeV. The results for the rank-zero vls 1 / 2 — > 
7r0p!/2 and v0d 3 / 2 ~ > 7r 0P3/2 transitions — the two al- 
lowed z^(1s0g?) — > 7r(0p) transitions — are given in Table 
H The (E x (j)) of Table H are calculated with the WBP 
interaction with the lowest allowed tiTiuj wave functions 
for each of the three nuclei involved in each of the four 
cases. In contrast to these simple wave functions, we 
present results in this study calculated for model spaces 
as complex as (l+3)7ia; — > (Q+2+A)fiui. The question 
might well be asked as to the relevance of Table [n] to 
these more complex — and hopefully more realistic - 
calculations. To address this question the procedure was 
applied to the vls 1 / 2 — » ffOpi/2 single-particle component 
in the 16 N — > 16 rank-zero decay with typical (l+3)7iu; 
and (0+2+4)?kj wave functions for 16 N(0") and 16 O(0+), 
respectively and with typical (0+2)?ia; wave function for 

the 15 N | states. In a (0+2)hcu space, there are 265 

15 N i states, 41 of which are spurious. The value of 
(E x (j)) found for the 224 non-spurious states is 32 keV, 
in close agreement with the value of zero keV associated 
with the simple calculation of Table O. 

For the ~> i"0f>3/2 transitions the large values of 

(E x (j)} mean that the matrix elements are not sensitive 
to its value, in fact, one might just as well (in ignorance) 
use HO wave functions for the j = § transitions and we 
do so for the calculation of the one-body matrix elements 
Mq and Mq and the two-body matrix element MJ and 

m;. 



S{p) = AE b (p)+E x -E f (8b) 

where Ei and Ef are the excitation energies of the initial 
and final states in the (A,Z— 1) and (A,Z) nuclei, respec- 
tively and AEb(n) is the difference in binding energies of 
the ground state of the initial nucleus and the (A— 1,Z— 1) 
core nucleus, i.e. the neutron separation energy for the 
ground states, and similarly for AE b {p). 

The main effect of small separation energies is to in- 
crease the relative magnitude of the "tail" of the radial 
wave function. As the separation energies increase, the 
difference between WS and HO wave functions lessens, 
and becomes insignificant compared to our knowledge 
of the wave functions. Thus we use HO wave functions 
for all single-particle transitions except those allowed for 



B. j3 and n rates 

Note that all the numerical results given in this subsec- 
tion are relevant to 16 N <-> 16 0. The relationship between 
experiment and theory is given here in a way that dis- 
plays the similarity between 0~ — > + (3 decay and + 
— * 0~ /i~ capture, it follows closely the treatments of 
Behrens and Biihring [jllj for decay and Nozawa, Ku- 
bodera, and Ohtsubo (NKO) (3lJ for both (3~ decay and 
fi~ capture but especially the latter: 

A/ 3 = |^K| 2 , (9a) 
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A„ = C, 



G 2 



' ^ R 2vl+u/M f A Ce 2 



9A(q 2 ) 
9A(0) 



\Mg\*. 



(9b) 



All quantities are in natural units h = c = m e = 1. The 
unit of time is the second and of length Xc e — the elec- 
tron Compton wavelength divided by 2ir. The decay rate 
A (= 1/t) has the units of sec -1 . G is the standard weak- 
interaction constant. The /3-decay phase-space factor is 



"o 



C(W)F(Z,W)(W' t - 1)5 W (Wo -WydW (10) 



where C(W) is the shape factor, F(Z,W) is the Fermi 
function, W is the electron energy and Wo the total dis- 
integration energy — both including the rest mass. The 
allowed phase-space factor, /o is given by the integral 
of Eq. @ with the shape factor C(W) = 1. For a 



pure rank-zero (R0) decay as is involved here, C(W) = 
f/fo = | Mo 3 1 2 where the M Q (a = /3,/i) of Eq. (|) are 
combinations of matrix elements. We choose to give all 
matrix elements in /to; thus Xce appears in Eq. (||). 
The axial-coupling constant <7a(? 2 ) is a function of the 
four-momentum transfer q. For nuclear j3 decay q = 
to a good approximation. We incorporate <7a(0) (= (ja) 
into the matrix elements. Hence it does not appear in 
Eq. (9a) and normalizes the g-dependent ratio gA(q 2 ) hi 
Eq. (|9b|). In Eq. (pb|), Cr is a second-order relativis- 
tic correction (see Appendix A), lu is the muon-neutrino 
energy, (1 + uj/Mf)^ 1 is a recoil correction with Mf be- 

2 

ing the mass of the final nucleus, and |0is(r)| r _ o is the 
probability of finding the at the origin (see Appendix 
A): 



|01 S (O)| = ft||<M0) I Point =H% 



(Zamlf 



(11) 



nucleus 



InEq. @ 



is the reduced \x mass and TZz is a 



correction factor obtained by solving the Dirac equation 
for the wave function of the muon in the field of a finite- 
charge distribution (see Appendix A). Parameter values 
(in natural units unless otherwise specified) relevant to 
16 N(0") 16 O(0+) are given in Table Using these 
parameters we find 



A/3 = 1.4206 • 10- 4 |M T = 0.485 ± 0.019 sec -1 , (12a) 



A, 



0.1210 \Mg\ 2 



1560 ± 94 sec" 1 



(12b) 



where the experimental results on the right are taken 
from Refs. || and |57|. The experimental rates lead to 
the following experimental matrix elements: 



= 83.90[A /3 (s- 1 )]2 = 58.4±1.1 fr 



MX 



2.874[A Al (s~ 1 )]2 = H3.5±3.4 /to. 



(13a) 



(13b) 



Results for the /3-decay matrix elements for the other 
three cases of interest are given in Table §. It is these 



experimental matrix elements that we will compare to 
theory. We write Mq and Mq in the form 

M* = [M? + M@ + a^M s } = [M?eP + a" s M s ] (14) 



[a%M T + M% — a^M s ]= [a!iM T e'" - a^M^' 



where the M" are the two-body meson-exchange cur- 
rent (mec) matrix elements defined in Sec. III.C below. 
When Eqs. (|l3|) and ( O ) are used to compare experi- 
ment and theory, it is conventional to extract values of 
e a that are required to reproduce experiment, and these 
will be referred to as e^ xp . The calculated values for e a 
based upon these equations to gethe r with the calculated 
two-body M" discussed in Sec III C| will be referred to as 



Our evaluation of the /3-decay rate follows the rigorous 
and accurate treatment of Behrens and Buhring |ll[| . In 
writing the expression for Mq in Eq. ( |l4| ) small terms 
included in the first-order treatment have been neglected. 
However, for the decays in question these terms con- 
tribute less than 0.2% to Mq 3 . In the Behrens-Biihring 
treatment the of Eq. (|lj) is given by |pT|,p~ 



38 



(Qp + 1) + = 7.208 + 3.187rg, 



(15) 



where Qp is dimensionless (in units of the electron mass) 
and the numerical results apply to the 16 N — > 16 decay. 
Note the r$ and the and Oj. of Eqs. (14) and (15) 
are positive-definite quantities so that the contributions 
of Mq and Mq add destructively in forming Mq and 
constructively in forming Mq. 

The a£ and a£ of Eq. (O) and r P s of Eq. @ are 
defined as ratios of matrix elements evaluated with extra 
radial factors to the normal matrix elements. The radial 
dependencies for /x~ capture are given explicitly in Ap- 
pendix A. In Eq. ( p5| ) the ratio is insensitive to the 
wave functions used. Approximate values of are f=a 
7.80, 8.52, 7.45 and 9.40 for the decays of n Be, 15 C, 16 C 
and 16 N respectively. Note however that all dj. and 
are calculated explicitly for the wave functions at hand. 



C. The two-body soft-pion contribution 

Towner j3j| has recently made an investigation of the 
one-pion exchange contribution to single-particle matrix 
elements of the axial charge in (closed shell ± 1) nu- 
clei from A ~ 16 to 208. Towner's results indicate that 
for A ^ 16 the soft-pion diagram alone gives an 

adequate representation of the meson exchange and so 
we only consider this term. The soft-pion contribution 
was incorporated into the shell-model calculations in the 
manner described for the similar parity-nonconserving 
(PNC) matrix element |42[ . The general expression for 
the soft-pion contribution is a sum of the product of a 
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(2) 

two-body transition density Dq \j1j2JTj3j4JT ) and a 
two-body meson-exchange matrix element: 



Mf = G xnffv G 



softiv^TJ 2_^i 



-D ( *\j lj2 JT] 3 3 4 jr) x 



(16) 



0\ <32 
J3<J4 
JTT' 



(iii 2 JT|||g(r r )i(o-i + <t 2 ) • r r (n X T 2 )Y^)\\\j 3ji JT') 

where x v — m^rr with r r = r± — r 2 , Y 7T (x 7T ) = (1 + 
l/x 7r )e" X7T /x^ and g{r r ) is a short-range correlation 
(SRC) function. The SRC used in this work was g(r r ) 
= 1 — jo(<7c r r) with q c = 3.93 /m _1 We will also 

compare results obtain with this SRC and with the sim- 
ple cut-off factor 9(r r — d) for which g(r r ) = 1 for r r > 
d and for r r < d. The soft-pion coupling constant is 
defined as 



V2gl 



NN 



8tt g A M 



2 A 

N 



Ce 



69.74 fm 



(17) 



where we have used g^NN = 13.684 |g2j. The two-body 
transition density of Eq. dlq) is given by 



D^(jij2JTj 3j4 JT') = 



(18) 



was made in Eq. (16) to obtain M% 



The relative contributions of M@ to Mq is formulated 
in terms of the ^ parameter defined in Eq. ( p^ ) . It is in- 
formative to consider the values of e^ ec obtained for sim- 
ple configurations. Evaluations with HO and WS wave 
functions are collected in Table |l^ for two different SRC 
factors. It is seen that ef„ ec has a fairly strong state de- 
pendence. Because of this state dependence and because 
the two-body and one-body matrix elements have differ- 
ent dependencies on the nuclear wave functions, Eq. jl^ ) 
or an equivalent expression must be used for a rigorous 
evaluation of e(^ ec . However, a useful approximation is 



E 3l3f d [ r ] i.y,y/H,:„,;./,y/).uJ (./,j/) 



(20) 



with a similar expression for e£ 



Comparison shows 



that this approximation — used with the results listed in 
Table IV to represent the ^medjijf) — underestimates 
e mec by ~5%. Nevertheless, this approximation is useful 
because the evaluation of Eq. (|l6|) takes an unusually 
large amount of computer time. In comparisons of dif- 
ferent diagonalizations and tests of different radial forms 
or SRC factors we make use of Eq. (EQ) . 



[(2AJ + 1)(2AT + 1)(1 + 5 jlh )(l + %,J]V2 • 

As for the other shell-model calculations, the evaluation 
of M@ involved four oscillator shells and the use of mixed 
HO and WS radial wave functions. Our results for single- 
particle transitions are identical to the soft-pion results 
of Towner [^9[ . Here we consider more complicated tran- 
sitions. The computer program used is formulated in 
terms of harmonic-oscillator wave functions. Our method 
of allowing Woods-Saxon radial wave functions is to ex- 
pand the appropriate WS radial wave function in terms 
of HO wave functions (up to 10 terms). As explained in 
Sec. III. A, we use WS wave functions for the vlsx/2 ~ ¥ 
7r0pi/2 transition but use HO wave functions to evaluate 
the rest. As for the one-body matrix element, cap- 
ture differs from 8 decay in that the analogous M% to 
the M@ of Eq. ( |l6| ) should be evaluated with the factor 
a T [~ Jo (<*"")] of Eq. (A17) of the Appendix inserted. 
Thus a dependence on r [= |(ri + r 2 )] as well as on 
ty is introduced. This dependence is approximated by 
retaining terms up to order to 2 with the result that the 
replacement 

(<Tl + <T 2 ) ' TV -> (CT! + <7 2 ) ' T r 

-l w 2 (r 2 + Ir^K + o-^-rv (19) 
+ \u 2 {r ■ r r )(cr 1 - <r 2 ) • t> 



IV. RESULTS FOR 16 N((T) <-> 16 O(0+) 

As a first orientation, consider the four decays of Table 
I in successive degrees of model-space complexity. The 
simplest of these is the "single-particle" approximation 
in which only the Opi/2 and lsi/ 2 degrees of freedom are 
allowed. Next allow for mixing within a single major- 
oscillator shell but keep only the lowest possible Tilo con- 
figurations, e.g., lhcu <-> Ohio for the 16 N <-» 16 O tran- 
sitions. In this approximation the WBP and WBN in- 
teractions are equivalent. The results for the two calcu- 
lations are compared in Table [v| First we observe that 
in all cases there is a significant reduction in going from 
the simple configuration to the major-oscillator configu- 
ration. This is primarily due a destructive interference 
between the v\s\ji <-> 7r0pi/2 an d v0d 3 /2 «-> ^Qpz/2 terms 
in Eq. (|. 

The next level of approximation allows for those terms 
which enter in the next order in a perturbation expan- 
sion. For lfko — > Ohio transitions, this means including 
2Uu> admixtures in the final state since these connect di- 
rectly to the Vtko initial state; Zfioj admixtures in the 
initial state are not included in this order because they 
do not directly connect to the Ohio final state and it is 
inconsistent to include them without also including ATiuj 
terms in the final state. A consistent higher-order calcu- 
lations involves (l+3)?iw — > (0+2+4)7iw transitions and 
is only possible at present for 16 N decay, which we now 
consider. 

The model-space dimensions for 16 N 0~ in a (l+3)7kj 
model space and 16 + in a (0+2+4)?kJ model space 
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are 713 and 4255, respectively. The gap reduction for the 
WBP interaction is 3.0 MeV. The lowest five 16 0+ T = 
states are calculated to lie at 0.00, 6.59, 10.77, 11.69 and 
14.25 MeV with the 6.59-MeV state predominately (88%) 
Ahui. The structure of the ground state is 43.2 percent 
OTuj, 43.2 percent 2hu> and 13.6 percent 4%ui. The lowest 
two 0" T = 1 states occur at 12.38 and 17.37 MeV and 
the structure of the lowest 0~ and 2~ T — 1 states are 
both (coincidentally) 71.1 percent lhu) and 28.9 percent 
3Tluj. 

A partial WBN level scheme is compared to experiment 
in Fig. ||. The gap reduction for the WBN interaction 
is 3.5 MeV. The structure of the ground state is 34.9 
percent Ofrco, 46.7 percent 2Hlu and 18.4 percent 4Kui. 
The 6.44-MeV 0^ state is predominantly (87%) Afiui. The 
structure of the lowest 0~ T = 1 state is 64.5 percent ITlui 
and 35.5 percent 3huj. 

The (3~ R0 and R2 decays of the 0~ and 2~ states of 
16 N are illustrated in Fig. |. Results for the 16 N(0~) -> 
16 O(0 + ) R0 (5 decay are given in Table [v| and those for 
the 16 O(0+) -> 16 N(0") R0 [T capture are given in Ta- 
ble VII. As a further test of the wave func tions, we com- 
pare with experiment Hffi} in Table VIII R2 results for 
the 16 N(2~) 16 O(0+) unique first-forbidden (3 decays, 
which are simpler and better understood (l^j2^] than the 
R0 decays. For example, no appreciable mesonic contri- 
bution to R2 transitions is expected. As in Ref. the 
good agreement between theory and experiment for these 
R2 transitions provides some evidence for the correctness 
of our wave functions. 

The decompositions of M^ and Mf obtained with the 
WBP interaction for the R0 and R2 f3 decays to the 16 O 
ground state into the four possible whuj — > (n±l)?ia; com- 
ponents are shown in Fig. ^. The various nhuj contribu- 
tions follow the classic pattern found for de-excitation of 
El-like particle-hole configurations in a previous study of 
first-forbidden decays in the A = 40 region f24| ]. In par- 
ticular, the Ihw — ► Ofku and 3Hlo — > 2fko matrix elements 
and the lfku — * 2Huj and 3huj — + Ahto matrix elements are 
closely equal and the latter two are out of phase with the 
former. The detailed composition of the Mq of Eq. (W) 
for the 16 N(0~) ->■ 16 (0+) f3~ decay is given in Table pa 



A. 15 C(| + ) 15 N(±~) 
The dimension D(J=i) of the J-matrix for the T 



T7T 1 

J — 2 



states of 15 C is 2369 in the four-shell 



2_ 

2' J 



states 



(l+3)?iu; model space. For the T 

of 15 N D(|) is 23762 for a (0+2+4)7iw calculation and 
265 for a (0+2)?ia; calculation. We can easily manage 
the 15 C (l+3)7ia> and 15 N (0+2)hw calculations but the 
(0+2+4)?ia; calculation is beyond our present resources. 
An attempt to truncate the 15 N (0+2+4)7ia; calculation 
by restricting the ATiuj part to 4p-4h (4 particle-4 hole) 
excitations between the Op and IsOd shells failed because 
the low-lying states in this truncation were highly spuri- 
ous. This is in contrast to similar calculations performed 
for the 0+ states of 16 |l3) and 12 C |lj. ( The low-lying 
+ states in the latter calculations were nearly spurious 
free.) However, the importance of the ATiuj component 
in the 16 O ground state suggests a similar importance in 
15 N and we would like to include an estimation of its ef- 
fect. Thus a (l+3)7iu; — > (0+2)?iu> calculation was made 
with the effects of adding a &huj term to 15 N estimated by 
assuming the same amplitude of Ahu> in the 15 N ground 
state as in 16 O (see Fig. 4) and assuming 

M(3hcj -> 4hcj) » M ( 1 ^^ 2 M M(3?kJ ^ 2huj) (21) 

M(Yhh> — > Ohuj) 

where M[(n ± l)^ — > iAlo] is the specific component, 
such as in Fig. 4, of either or M§. Thus Eq. <M) 
quantifies our observations on the systematics of Fig. 4. 

A value of A^nuj — —8.50 MeV was determined for the 
15 C (l+3)hoj calculation by a least-squares matching to 
the eleven experimental even-parity A = 15 T = | energy 
levels for J < | states below 11 MeV, the yrast T = § 



£ , and § states, and the yrast T 



1 13 + 

2 2 



state. The 



value of A 2 htu used in the 15 N calculation was determined 
as —8.30 MeV from a similar consideration of the A = 
15 odd-parity spectrum. 

The results are Mj = 27.66 /to, Mf = 1.53 /to, 4 xp 
= 1.54 ± 0.04 (uncertainty from experiment only). 



V. RESULTS FOR 15 C, "Be, AND 16 C /T DECAY 



In this section we will only consider results obtained 
with the WBN interaction using the gap method and 
mixed WS and HO single-particle radial wave functions. 
We will confine our comparison of experiment and the- 
ory to a quotation of the Mq and the value of e^ xp , which 
when combined with the calculated Mq reproduces ex- 
periment. In all three cases, the calculated e^ ec differ 
negligible from that calculated for 16 N decay. 



1 Be(i + ) 



l B(|") 



The calculation for llBe( 



11B( 



follows the 



procedure just described for 15C(^ ) — ► 15N(i ) with 
the reliance on the (0+2+4)Sa; 16 Of wave function 
in the latter case changed to a similar reliance on the 
(0+2+A)hio 12 C Of wave function. The D(±) for the 
(l+3)7iu; — > (0+2)huj calculation are 5674 and 1063, re- 
spectively. Little of a definite nature is known about the 
2hu> and 37iw states of n B and 11 Be. Thus the values of 
A.2huj and A^h^ used in the calculations were assumed to 
be equal and were obtained by linear interpolation be- 
tween values found for 10 B, 12 C and 15 N. The resulting 
value is -4.00 MeV. 
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The effect of the ATilo component in n B was consid- 
erably less than the similar effect in 15 N simply because 
the 12 C ground state wave function has only 1.86% ATiu) 
as opposed to 13.6% in 16 0. 

The results are = 14.05 fm, Mjj = 1.10 fm, e^ xp 
— 1.56 ± 0.08 (uncertainty from experiment only). 



experiment, e" . It is evident from Tables VI, and VII 



5 C(0+) 



3 N((T) 



After considering 16 C(0 + ) — > 16 N(0~) in detail we con- 
clude that the calculated matrix elements are unusually 
sensitive to the details of the calculation and thus not 
suitable for a determination of e@ xp . To understand this, 
first consider the simple 2huj — > ITiu calculation of Table 
[v|. In the case of 16 C the reduction of the Mq in go- 
ing from the simple single-particle configuration to the 
major-oscillator configuration is unusually large and in- 
teraction dependent, because of the competition between 
the 0d^ 2 and ls 2 / 2 configurations in the 16 C ground 

state. As a consequence of this, the configuration in 16 C, 
which is responsible for the destructive vOdz/2 — > tt^Pz/2 
term (but which accounts for only about 0.3 percent of 
the 16 C wave function), results in a ~50% reduction of 
the matrix elements. 

Now consider the mixing of 2Hlo and Afku components. 
A (2+A)huj calculation — which has D(0+) = 4055 — was 
performed. The same decrease of the Op-lsOd energy gap, 
3.5 MeV, for the T = 2 + states was used as was applied 
for the T = + states. We find the first 10 states (i.e., 
all that were examined) to be >44% Ahcu with, e.g., the 
ground state being 55% 2Huj and 45% 47kt>. The (0^,2) 
excitation energy in 16 was 22.2 MeV as compared to 
the experimental value of 24.52 MeV. The prominance of 
the Ahoj component may seem surprising at first sight; 
however, if a weak-coupling approximation such as that 
of Bansal and French J29| is used to estimate the binding 
energies of 14 C® 18 0, and 12 C<g> 20 O, one finds that these 
energies are nearly degenerate so that which lies lowest 
and the energy gap between them depends on the details 
of the interaction. Thus it is not surprising that the 
{2+4)Tiuj mixing is large and unusually sensitive to the 
details of the calculation. 

Because of these sensitivities we believe that the 16 C 
decay is primarily a test of the wave functions and does 
not provide a good measure of the mesonic-exchange- 
current enhancement. This extreme wave function sen- 
sitivity for 16 C also applies to similar R0 decays of 17 N 
and 18 Ne We do not consider this decay any further 
in this study. 



VI. DISCUSSION AND SUMMARY 

We comment first on the value of the enhancement fac- 
tor for time-like axial-charge matrix elements deduced 
from a comparison of the impulse approximation with 



that the value of e" xp is strongly dependent on the single- 
particle radial wave function used to describe the vOsi/2 
— > transition. The WS result is strongly preferred 

since it provides the most realistic estimate of these ra- 
dial wave functions. The HO results are included to give 
some indication of the sensitivity of the results to the ra- 
dial form. We will not consider the HO results further. 
Likewise, the results of the WBN interaction are strongly 
preferred over those of the WBP interaction for the rea- 
sons stated in Sec. II. Again, the WBP results are listed 
to give an indication of the sensitivity to the shell-model 
interaction and will not be considered further. 

In addition to the dependence on single-particle wave 
functions and shell-model interactions, there is a further 
dependence we have not yet considered; namely, the de- 
pendence of the fi~ capture rate on the pseudoscalar cou- 
pling constant g p (see Appendices A and B). For the 
WBN interaction with WS wave functions the depen- 
dence on g p can be expressed as 



1.55 + 0.085(p p - 6.939). 



(22) 



In the evaluation of the (3~ R0 matrix element Mq the 
contributions of Mq and Mq add destructively while 
they add constructively in Mq [see Eq. (Q)]. This rela- 
tive phasing results in less model-dependence for the de- 
termination of e^ xp than for that of e@ xp and it has often 
been asserted that, for this reason, /i~ capture provides 
a considerably more reliable determination of e exp than 
(3 decay. However, the uncertainty in g p (see Appendix 



B) completely negates this conclusion; from Eq. 
see that the range of g p of ^7-12 allowed by analysis of 
radiative fi~ capture (see Appendix B) corresponds to a 
range in e exp values of 1.55-2.00 — a range large enough 
to offset the afore-mentioned advantage. 

Our adopted value for e exp follows from a consideration 
of the results for the three (3~ decays that were analyzed. 
The average of the three results (1.54 ± 0.04, 1.56 ± 0.08, 
1.63 ± 0.02) quoted for e exp is 1.60 ±0.03 and this is the 
result we shall adopt. We note that the fi~ capture result 
of 1.55 ± 0.08 (for the PCAC value of g p ) is in good 
accord with this value especially if the small difference 
between and e£ of ^0.02 is added to the e^ xp value. 

The previous most ambitious calculation of the pro- 
cesses in question were the (1+3)Tiuj — > (0+2)fiLu calcu- 
lation of Warburton |45|] for the four f3~ decays of Table 
I and the (l+3)huj <-> (0+2+4)?ki; results of Haxton and 
Johnson Q for the 16 N(0") <-> 16 O(0+)/3- and fj,- pro- 
cesses. The calculations of Warburton [[l5|, using the MK 
interaction, result in an average e@ xp of 1.64 for the de- 
cays of Table I, in good agreement with the present result 
of 1.61. Haxton and Johnson obtained results very close 
to the present results given for the WBP interaction with 
WS wave functions . This is not surprising to us since 
we had found that the V 2 ^ derived from the Kuo bare 
G matrix Q — the V 2huJ of the MK interaction used 
by Haxton and Johnson — has quite similar properties 
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to the V 2huJ of the WBP potential. It is the change to 
the Bonn potential for V 2huj that gives us a substantial 
improvement in the agreement of experiment and theory. 

The theoretical values for the enhancement factors, 
e mec an d e meci calculated from a meson-exchange model 
in the soft-pion approximation are very close to our 
adopted experimental value and show very little model 
dependence varying from 1.60 to 1.63 and 1.59 to 1.61, 
respectively, depending on the choice of interaction and 
single-particle wave functions. This insensitivity has 
been noted previously in work on parity non-conserving 
interactions and follows because the two-body ma- 
trix element can be represented quite well by an effective 
one-body matrix element proportional to a-p — and thus 
to Mq — consequently the ratio of the two-body to one- 
body matrix elements is not sensitive to nuclear struc- 
ture. However, as discussed by Towner [^9), the results 
from meson-exchange models are dependent on the choice 
of a short-range correlation function. This is particularly 
true for the short-range operators discussed in |}9| origi- 
nating in heavy-meson exchange. In this paper, we have 
only considered the long-range pion-exchange operators 
evaluated in the soft-pion approximation for which the 
sensitivity to SRC is somewhat less. Some sample calcu- 
lations are given in Table IV. For the all- important vlsi/2 
— » irOpi/2 transition the use of the 9(r r — 0.71) form re- 
duces e^ ec by a factor 0.96 compared to the 1— jo(3.93r r ) 
form and the other transitions in this table have similar 
ratios. In the full calculation for both /3-decay and \i~ 
capture (a = (3,fi), the reduction factor is 0.95-0.96. In 
what follows we will only discuss calculations using the 
SRC function g(r r ) = 1 — jo(q c r r ) with q c = 3.93/m -1 
II- 

As discussed above and in Appendix B, the jjT capture 
results depend on the pseudoscalar coupling constant, g p . 
This dependence offers the opportunity to extract a value 
for g p from a comparison of the (3~ and /i - results. For 
the WBN interaction and WS wave functions with the 
l—jo(q c r r ) SRC — our preferred selection — the adopted 
value = 1.61 ± 0.03 (which we assume to yield e£ 
= 1.59 ± 0.04) and use of Eq. (22) leads to g p = 7.4 
± 0.5 consistent with the PCAC value. Note that the 
uncertainties assigned to e exp and g p reflect experiment 
only and do not include any estimate of the uncertainty 
associated with the theoretical analysis. 

Finally we compare these results for the A = 16 mass 
region with our recent results for the A = 132 and A 
= 208 H regions. As noted in Sec. III. A the calcula- 
tions in the heavier mass regions arc further complicated 
by the need to introduce core-polarization corrections. 
These corrections depend sensitively on the strength of 
the residual interaction used in their evaluation and in 
particular on the strength of the tensor force which, ac- 
cording to Brown and Rho p8| , is liable to medium mod- 
ifications similar to those proposed for the pion-decay 
constant fljjffiOj . With a weak tensor force, the enhance- 
ment factors deduced are e% x (A=132) = 1.82 ± 0.07 



and ef a;p (A=208) = 1.79 ± 0.04. With strong tensor 
forces, values of e^ xp some 10% larger are obtained. Thus 
there is clearly more enhancement in heavy nuclei than 
in light nuclei, which are characterized by our current 
value e exp (A=16) = 1.61 ± 0.03. To quantify the mass 
dependence we define a ratio, 

_ e exp (A = 208) - 1 
r e exp (A=W)-l> [Z6) 

which for weak tensor forces has a value 1.30 ± 0.10 and 
for strong tensor forces a value 1.44 ± 0.11. The cal- 
culated value of this ratio from meson-exchange models 
obtained by Towner p9[ is 1.38, which is approximately 
midway between the experimental results obtained as- 
suming weak and strong tensor forces. There is still some 
room for further sources of enhancement in the heavier 
nuclei such as medium modification of the nucleon mass 
and of the pion-decay constant used in the soft-pion ap- 
proximation as proposed by Kubodera and Rho |5(| on 
the basis of a scaling in the effective chiral Lagrangian 
discussed by Brown and Rho p9| |. But the compelling 
need for such an alternative explanation that was thought 
present a few years ago now appears to be absent. 
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APPENDIX A: yT CAPTURE 

Nozawa, Kubodera, and Ohtsubo |3l| — referred to as 
NKO — developed an approach to fi~ capture on le O, 
which can be used with numerical solutions of the Dirac 
equation. In this appendix we wish to make the bridge 
between the /i~ capture formalism developed by NKO 
and numerical results found by combining their solution 
of the Dirac equation with our shell-model wave func- 
tions. In essence this means we will give explicit expres- 
sions for the radial dependence of the matrix elements, 
which will then be evaluated by numerical integration. 

Our expression for the [i~ capture rate [Eq. (j9~b|)] is 
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a -r G2 ^ ^ ls(0)|; 



gA(g 2 ) 

Sa(0) 



-. 2 



l M o I 2 - 



(Al) 



In order to evaluate relativistic effects, NKO applied a 
Foldy-Wouthuysen transformation first to order 1/Mjv 



with the coefficients given in Table [|. 

The Foldy-Wouthuysen transformation to first order in 
I /Mm leads to a one-body operator for \i~ capture given 

by m 



and then to order l/M 2 N . They showed that the terms ^ = t* A ' « + M + + (A8) 



AT 

of order 1/Mjy contribute ~ 2% (constructively) to the 
H~ capture matrix element Mq. We will calculate A M to 
first order in 1/-Mjv b ut i nclude a multiplicative factor 
of Cr (=1-04) in Eq. flAl] ) to compensate for the 1/Mjy 



term in our evaluation of Mq . 

In Eq. jAlj ) |^i7(0)| 2 is tne probability of finding a K- 
shell muon at the origin and is just 1/4tt times the square 
of the large component of the muon wave function G_i(r) 
evaluated at r = 0. For a point nucleus flit] 



G9(r) = 



(1+7)A 



r(l + 2 7 ) 



1/2 



L (2Ar)' 1 



(A2a) 



with radial functions 



L±(r) = -=[G_i(r)ji(wr) ± ^(r^wr)], (A9) 



V2 



1 



£±(r) = [G_ 1 ( r )io(o;r) ± F_i(r)ji(u/r)], (A10) 



V2 

and operators 



(All) 



a2 
1 + 7 



(A2b) 



where -F_i(r) is the smaller component of the muon wave 
function and 7 = [1 — (aZ) 2 ] 1 / 2 . When appropriate, 
we use the superscripts (p) and (gr) on G_i and to 
denote a point nucleus and a nucleus with a Gaussian 
charge distribution, respectively. In Eq. (A2) 



A = 



K + £ M )(m;-i=y 



1/2 



= mil (A3) 



where is the energy eigenvalue for the K-shcll for a 
point nucleus. The approximation often used for (r) 



G W( r ) ^2{aZmlf' 2 . 



(A4) 



Thus we take the probability of finding the muon at the 
origin (in the unit volume A Ce ) to be 



l<M0U 



4tt 



-(aZm 



r \3 n2 



(A5) 



where Rz is a correction evaluated by solving the Dirac 
equation for a realistic extended nuclear-charge distribu- 
tion. 

We use NKO's solution of the Dirac equation. These 
authors used a Gaussian charge distribution p(r) given 
by 



P( r ) = Po[l + a(r/r o y 



-(r/ro) 2 



(A6) 



with ro = 1.83 fm and a — 1.45. This charge distribution 
gives G_l(0) = 0.91688 and large and small components 
of the muon wave function, which are well reproduced by 
the phenomenological power expansion 
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G^(r) = 5>„, 



n=0 



10 

71 = 



F^(r) 



b n r n , (A7) 



A = ig A a ■ P/2M, 



N, 



A = m^GpCr ■ k/2M 



N- 



(A12) 



(A13) 



Here ja(ojr) and j\ (tor) are spherical Bessel functions rep- 
resenting the neutrino wave functions. Note that P and 
k depend on the momenta of the initial and final nu- 
cleons, P = pi+ pf and k = pi — p/, which occur in 
combination with radial functions C(r). In the Fourier 
transform to coordinate space these momenta transform 
into derivative operators according to the replacement 
rules: 



P£(r) -» -iiJLcM - 2i£(r)V, 
Or 



k£(r) — * ir— £(r). 
or 



(A14) 



(A15) 



Then, upon evaluating the partial derivatives, we find 
expressions for the full matrix element evaluated to order 
l/Mjy. In terms of the operators defined in Eq. (1) and 
the normalization imposed by the expression for the /i _ 
capture rate in Eq. (Al), the one-body operator becomes 



(A16) 



with 



1 



G-i(0) L 



G_i(r)j (wr)-.F_i(r)ji(wr) , (A17) 



1 



rG_i(0) 



[giarge + 8i arge (r)]G -i(r)ji(ujr) (A18) 



\i3 small H~~ $small (r)]F_i(r)jo(ur) 
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where in our notation Mq and Mq are equivalent to 
NKO's notation — -ri— a ■ V r_ and a ■ r r_. In Eq. 
(A18) 

fflarge = -g» = 1 - (&> - 1) = 0.6969, 



= 1 



2M W 



^large 



(r) 



2m„(<? p - 1) + cj(. 9p + 1) = 2.0753, 



-[y(r)+^]( 5p + l) 



2M 



in 



5 smaller) = 



2M 



N 



[V(r)+E K }(g p -l) + -(g p + l) 



and g p = m )1 G p {q 2 ) j 'g A {q 2 ) is the pseudoscalar coupling 
constant. In these expressions, Ejc is the iVshell fi~ 
binding energy [= (1 — 7)m M for a point nucleus], and 
V(r) is a potential evaluated from the charge distribution 
used to calculate G_i(r) and F_i(r). The potential is 



V(r) 



1 



p(x)x 2 dx + / p{x)xdx (A19) 



with 



p(r)dr = a2. (A20) 



The approximate solution designated "simple" by 
NKO corresponds to the limits 



G_i(r) -> GLi(O), 
F-i(r)->0, 
F(r) + - 0. 



"simple" solution (A21) 



We have arranged the terms in Eqs. (A17,A18) so that 
it is easily apparent that in this limit Eq. (A16) reduces 
to 

Jl- bod y (simple) = j (wr)M T + 9ll —ji (wr)M$ . (A22) 

P LOT 

The terms in & small (?*) and 6i arge (r) give small contribu- 
tions (< 2%). As pointed out by NKO, the g sma ii term 
gives a contribution of order 10%. 



APPENDIX B: PSEUDOSCALAR COUPLING 
CONSTANT 

For semi-leptonic weak processes, the V — A structure 
of the weak interaction is modified by the induced weak 



currents that arise from the presence of the strong inter- 
action. For the axial current one common parameteriza- 
tion is 



= iu(pf)[g A (q 2 )7^75 



G P (g 2 
2M 



Qp, il^uiPi) 1 ^ (Bl) 



where the strong interactions have renormalized the 
axial- vector coupling constant gA (<? 2 ) from a bare value 
of unity and introduced a pseudoscalar term with cou- 
pling constant G p (q 2 ). Here u(j>f) and u(j>i) are spinors 
representing the final and initial state nucleons, q is the 
momentum transfer, and r a is the Pauli isospin matrix 
with the superscript a representing the Cartesian index 
x ± iy. If, additionally, a constraint is imposed on Eq. 
(Bl), namely that the axial current shall be partially 
conserved, the PCAC condition, then a relation between 
G p (q 2 ) and g A {q 2 ) can be established. 

One model for the weak axial current of a nuclcon, 
that satisfies by construction the PCAC condition, is 
one of meson dominance discussed by Towner [^9|. In 
this model, the g A term in Eq. (Bl) originates in the 
axial current being mediated by the A\ meson in its 
interaction with a nucleon, while the G p term is medi- 
ated by the 7r-meson. The chiral Lagrangian that is used 
to describe the meson-nucleon interaction preferentially 
chooses pseudovector coupling for 7rNN vertices. Hence 
there is a pseudovector form, g M ^75, to the G p term in 
Eq. (Bl). However for on- mass-shell nucleons, the use 
of the Dirac equation can transform this equation into 
one with pseudoscalar coupling, ^75. Using this meson 
dominance model, Towner p9|] obtains 



G P (q 2 



2Mg A (q 2 ) 



(B2) 



The last factor, which tends to unity in the limit 
/m A — > 0, is not usually present in the literature. To 
get this result from the chiral Lagrangian the Weinberg 
relation is used relating the Ai-meson mass to the p- 
meson mass: m 2 A — 2m 2 . With an experimental p-meson 
mass of 768.1 MeV, this leads to an Ai-meson mass of 
m A = 1086.3 MeV. On the other hand, neutrino-nucleon 
scattering analysed using g A (q 2 ) = g A (Q)m 2 A / (m 2 A + q 2 ) 
deduce that m A ~ 950 MeV. Thus there is some small 
uncertainty on the precise value of m A to use in Eq. (B2). 
We will take 771,4 = 950 MeV and use for the nucleon and 
pion masses their isospin averages to obtain for the di- 
mensionless quantity 



m p G p (q 2 
9A(q 2 ) 



2Mm p 
q 2 + m 



1 



6.939 (B3) 



0.8m,,. This is the value recorded in Table III. 



for q 2 

The pseudoscalar coupling constant, g p , is the least 
well measured of all the nucleon weak-interaction cou- 
pling constants. Its value in Eq. (B3) typically repre- 
sents a free-nucleon value and could be modified when 
that nucleon is embedded in a nuclear medium. The most 
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promising way to determine g p experimentally is in radia- 
tive [i~ capture f5lj . The branching ratio for radiative \x~ 
capture relative to ordinary (non-radiative) \jT capture 
is particularly sensitive to g p . However the extraction 
of g p from nuclear radiative \i~ capture measurements 
requires a model calculation of the inclusive nuclear re- 
sponse function and this piece of the determination is 
not yet under good control. For example the most pre- 
cise measurement for the le O branching ratio |5l[] yields 
values of g p = 7.3 ± 0.9, 9.1 ± 0.9, and 13.6 ± 1.9 when 
compared with three different recent calculations of the 
nuclear response. Nevertheless there is a hint here that 
the value of g p in le O is larger than the PCAC value, 
Eq. (B3). A similar statement can be made for 12 C. For 
heavier nuclei, analysed using a Fermi-gas model calcula- 
tion of the nuclear response, the indication ]5l] , |52]| is that 
g p falls below the PCAC value and even quenches to zero 
for nuclei as heavy as Pb. 

For le O, there is an alternative way to obtain g p ; 
namely, from the first-forbidden 0~ to + (3 transition 
from 16 N and the inverse /i~ capture on le O, the reac- 
tions under study in this paper. Gagliardi et al Busing 
the model calculations of Towner and Khanna [ |34| de- 
duce g p = 11 ± 2. This result is model dependent so that 
our evaluation — given in Sec. VI results in the quite 
different g p — 7.4 ± 0.5. Because of the model depen- 
dence involved, we have decided in this paper to quote 
results using the PCAC value of g p , but give sufficient 
details — such as Eq. (|2^) — that the results can easily 
be modified for different values of g p . 
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FIG. 2. Comparison of partial WBN and experimen- 
tal level schemes of 16 O. The experimental energies and 
spin-parity assignments are from Ref. j53|. The theoretical 
levels were calculated with the indicated downward shifts to 
the different nhu> components relative to the unshifted Qhoj 
components. The centroid of the lhu> + 3hu) shifts was set 
so as to minimize the difference in excitation energies for the 
yrast T = 1 odd-parity quartet. The other shifts were set as 
discribed in the text. The levels included in the figure are 
the lowest five + T = level, the odd-parity T = 1 quar- 
tet, and all T = odd-parity levels below 12-MeV excitation. 
There is no WBN counterpart for the known predominantly 
5hu) experimental 1~ and 3~ states at 9585 and 11600 keV, 
respectively. The experimental + level at 11260 keV is ten- 
tative. 153]. We assume it does not exist. 



FIG. 3. Experimental data for R0 and R2 (3~ decays of the 
- and 2~ states of 16 N. Each level is labeled by and E x 
(in keV). The R0 decay is labeled by the decay rate and the 
R2 decays is labeled by the total half-life and the individual 
branching ratios. 

FIG. 4. Schematic showing the contributions of nTiuj — > 
(n±l)ftw transitions to the 16 N <-> 16 R0 Mo matrix ele- 
ment in the 0~ <-> Of transition (4a) and the R2 M% matrix 
element in the 2~ «-» 0^ transition (4b). Both are calculated 
with WS wave functions and the WBP interaction. cti n and 
ctf n are the amplitudes of the various hui components in the 
initial and final states, respectively. 



FIG. 1. The four A J = decays of interest to the present 
study. Each level is labeled by J w and E x (in keV). Each /3~ 
branch is labeled by the percentage branching ratio (or decay 
rate) and log/oi value. 
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TABLE I. Experimental data for AJ = first-forbidden decays in the A ~ 16 region. ii/2_is 
the total half-life and b r the branching ratio for the indicated final state. Q/3 is defined in Eq. 
C(W) is the average shape factor, = f/fo [see Eq. (|l(J)], and is given by C(W) = Bj where 
B^' = |M^| 2 is the rank R (3 moment. For the | + — > | transitions R = 0,1 and the measured 
/B[ ^ listed in the sixth column is used to obtain Mq from C(W). For the + <-> 0~ transitions 
this step is not necessary. 



Transition 



s l/2 



6 r 
(%) 



(kcV) (/to) 



Ref. 



15 C(I + ) 



Bc(i + ) -» "B(|-) 

5 N(i-) 
16 C(0+) -> 16 N(0~) 
16 N(0~) -> 16 O(0+) 



13.81(8) 31.4(18) 9381.3(60) 14.4(4) 

2.449(5) 36.8(8) 9771.68(80) 32.8(4) 

0.747(8) 0.68(10) 7891.7(43) 13.4(9) 

1.429(56) 100.0(0) 10539.5(23) 58.4(11) 



<0.29 13.3(11) 

0.185(13) 29.6(11) 

0.00 13.4(9) 

0.00 58.4(11) 



TABLE II. Summary of results for the neutron and proton separation energies S(n) and S(p) for use with Woods-Saxon 
wave functions. All energies in the last seven columns are in keV. 



Transition 


j 


AE b (n) 


AE„(p) 


Ei 






S(n) 


S(P) 


16 N(0~) -» 16 O(0+) 


i 

2 


2491 


12128 


120 








2371 


12128 




3 
2 










6324 


8695 


18452 


16 C(0+) -» 16 N(0~) 


1 

2 


4251 


11480 





120 


25 


4276 


11385 




3 
2 










4267 


8518 


15627 


15 C(i + ) -> 15 N(±~) 


1 

2 


1218 


10208 











1218 


10208 




3 
2 










8746 


9964 


18954 


11 Bc(i + )^ 11 B(I-) 


1 
2 


504 


11228 





2125 


70 


574 


9173 




3 
2 










7013 


7517 


16116 



14 



TABLE III. Quantities needed in the evaluation of the 
16 N(0~) <->■ 16 O(0 + ) f3~ and /x~ matrix elements and decay 
rates. 



Quantity 



Definition 



Value 



fine structure constant 1 / 137.036 

electron Compton wavelength 386.159 fm 

uniform charge radius 3.537 fm 

aZ/2r u 3.187 

[1 - (aZ) 2 ] 1/2 0.99829 

Tz-avcragcd nuclcon mass 1837.41 

Tz-averaged pion mass 270.128 

muon mass 206.768 

reduced muon mass 0.99296m^ 
Weak-interaction factor (sec -1 ) [8851(20)] _1 

m - 186.144 

+ 2m M w 0.80051m 2 

four-momentum transfer 184.998 

Ga/G [= g A (0)] 1.261(4) 

[1 + (q/MN) 2 ]- 1 0.9900 

/3-decay Q value 20.624(5) 

phase-space .factor 1.2574 

See Eq. (U) 548.10 

See Eq. (11) 0.91688 

PC AC value 6.939 

7m^ 204.961 

(l-7)mj 0.350 

recoil correction 0.9937 

rclativistic correction 1.04 



Ac, 

7 
A/jv 



G 2 /2m 



<r 

Q 
9A 

9A(q 2 )/gA{0) 
Q/S 

TZz 

g P (q 2 )/gA(q 2 ) 

(point nucleus) 
E k (point nucleus) 
(1 + u/Mf)- 1 
Cr 



E. 



TABLE IV. Soft-pion enhancement factors for simple transitions encountered in 16 N(0~) — > 
16 O(0 + ) calculated with HO wave functions (Tuo = 13.60 MeV) for two different short-range corre- 
lation functions. WS results are also given for the vlsi/2 — * an d ^0^3/2 — * 7i~0p3/2 transitions. 
Note that the last three transitions have an 16 O Ohu core. 



Initial State 




Final State 


1 -io(3.93r r ) 


9(r> - 0.71) 


<4ec 

1 -jo(3.93t>) 


9(ry - 0.71) 








lhaj — > Ohuj 








t>ls 1/2 7r(0pi/ 2 ) 


-1 


(0s) 4 (0p) 12 


1.549" 


1.488 a 


1.510 a 


1.425 a 








1.523 b 


1.467 b 


1.509 b 


1.427 b 


v0d 3/2 ir(0p 3/2 ) 


-1 


(0s) 4 (0p) 12 


1.418" 


1.383" 


1.449 a 


1.387" 








1.428 b 


1.379 b 


1.457 b 


1.393 b 








ITiuj — •> 2huj 








vls 1/2 Tr(0p 1/2 ) 


-1 


f (lsi/aMOai/a) -1 


1.405 


1.367 


1.417 


1.361 


vls 1/2 Tr{0p 1/2 ) 


-1 


v(ls 1/2 ) 2 ir(0 Pl/2 )- 2 


1.525 


1.458 


1.488 


1.406 


iy0d 3/ 27r(0p3/ 2 ) 


-1 


^(0d 3/2 ) 2 7r(0p3 /2 )- 2 


1.331 


1.293 


1.356 


1.306 








IsOd -> Oflp 








0d 5/2 




O/5/2 


1.226 


1.201 


1.253 


1.212 


0d 3 / 2 




l?3/2 


1.488 


1.430 


1.436 


1.360 


lSl/2 




lPl/2 


1.305 


1.264 


1.310 


1.248 



a HO results. 
b WS results. 
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TABLE V. Results obtained for restricted model spaces. 
All were obtained with the Woods-Saxon potential. 



Transition 




Model space 


M 


M 


U Be(§ + )- 




single-particle 


34.3 


-1.59 




lhcj — > ohuj 


18.0 


-0.97 


15 c(i+)- 


15 N(i-) 


single-particle 


47.3 


-2.20 






lhcj — > Ohuj 


35.9 


-1.76 


16 C(0+) -» 


16 N(0~) 


single-particle 


81.1 


-3.12 






2huj — > lhaj 


26.8 


-1.22 


16 N(0~) -. 


16 O(0+) 


single-particle 


81.1 


-3.12 






lhu! — > Ohui 


62.1 


-2.70 



TABLE VI. Results obtained from the (l+3)hu; -» (0+2+4) Two calculations for 16 N(CT) -> 
16 O(0 + ) R0 P decay. The ef xp were obtained from Eq. ( |14[ ) using the experimental value of Mq = 
58.4 ± 1.1. The e^ec were obtained using the calculated matrix element as discussed in Sec. 
IILC with the [1— jo(3.93r r )] short-range correlation. 



Interaction WS/HO aj M| eg lee 

WBP HO 49~9 Tool -2.59 1.69(2) TTBT 

WBN HO 56.6 9.49 -2.26 1.41(2) 1.60 

WBP WS 41.7 9.37 -2.58 1.98(3) 1.61 

WBN WS 48.8 9.36 -2.26 1.63(2) 1.62 



TABLE VII. Results obtained from the (1+3)Tilu ^_(0+2+4)ftw calculations for 16 O(0+) -» 
16 N(0~) fj,~ capture. The e^ xp were obtained from Eq. ( J14| ) using the experimental value of Mq = 
113.5 ± 3.4. The e^ec was obtained using the calculated matrix element M£ as discussed in Sec. 
IILC with the [1— jio(3.93r r )] short-range correlation and the PCAC value of g p of Table III. 

Interaction WS/HO a^ ag M| e^" 

WBP HO WE 65~9 1.02(7) T760~ 

WBN HO 51.5 57.2 1.09(7) 1.59 

WBP WS 36.8 56.9 1.54(9) 1.61 

WBN WS 41.8 48.8 1.55(8) 1.59 



TABLE VIII. Results obtained from the (1+3)7Uj -> 



(0+2+4)?k<; calculations 


for 16 N(2") 


- 16 O(0+) 


R2 p decay. 


n Exp 


Interaction 


WS/HO 


Mf(0+) 


1 3.04(2) 


WBP 


HO 


2.91 




WBN 


HO 


2.61 




WBP 


WS 


2.93 




WBN 


WS 


2.63 


2 1.09(18) 


WBP 


HO 


0.96 




WBN 


HO 


1.03 




WBP 


WS 


0.96 




WBN 


WS 


1.03 



1G 



TABLE IX. WBN 16 N((T) -> 16 O(0+) results for the D^\jij s ) and matrix elements of Eq. (|) 
calculated with WS wave functions for the dominant line five and HO wave functions for the other 
entries as discussed in the text. 



vji 


tj7 












Op 1/2 


Os 1/2 


0.0019 


3.6473 


0.0069 


-105.7247 


-0.1988 


lp 1/2 


Os 1/2 


-0.0016 


0.0000 


0.0000 


-0.0187 


-0.0000 


Od 3/2 


Op 3/2 


0.0504 


6.9739 


0.3515 


-184.9045 


-9.3192 


Os 1/2 


Op 1/2 


0.0211 


3.6139 


0.0762 


105.8366 


2.2310 


Is 1/2 


Op 1/2 


0.7857 


-3.1164 


-2.4486 


73.0472 


57.3932 


Of 5/2 


Od 5/2 


-0.0062 


10.1179 


-0.0627 


-269.3061 


1.6697 


Op 3/2 


Od 3/2 


-0.0424 


6.9820 


-0.2963 


185.8391 


-7.8870 


lp 3/2 


Od 3/2 


0.0037 


-4.4158 


-0.0162 


117.5352 


0.4302 


Op 1/2 


Is 1/2 


-0.0243 


-3.1224 


0.0759 


-83.1269 


2.0216 


lp 1/2 


Is 1/2 


-0.0036 


4.9370 


-0.0178 


-131.4311 


0.4732 


Od 5/2 


Of 5/2 


-0.0003 


10.1179 


-0.0030 


269.3061 


-0.0808 


Od 3/2 


lp 3/2 


0.0013 


-4.4158 


-0.0057 


-117.5352 


-0.1516 


Os 1/2 


lp 1/2 


0.0000 


0.0000 


0.0000 


0.0187 


-0.0000 


Is 1/2 


lp 1/2 


0.0169 


4.9370 


0.0834 


131.4311 


2.2199 






Total 




-2.2564 




48.8014 



TABLE X. The power-series coefficients of Eq. (|A7|) that 
reprod uce t he Dirac wave function for the charge distribution 
of Eq. ( A6) with an absolute error < 5 x 10 -6 for r < 10 fm. 
The numbers in square brackets are powers of 10. 



n 


a>n 




b„ 







+ 9.16882[ 


-01] 


+ 0.00000]- 


-00] 


1 


+ 3.56365[ 


-06] 


- 8.29499[- 


-03] 


2 


- 4.51773[ 


-03] 


+ 3.49933[- 


-04] 


3 


+ 1.06662[ 


-04] 


- 3.02963[- 


-04] 


4 


- 5.80056[ 


-05] 


+ 3.88561[- 


-04] 


5 


+ 7.30971[ 


-05] 


- 1.40780[- 


-04] 


6 


- 2.20053[ 


-05] 


+ 2.59785[- 


-05] 


7 


+ 3.29970[ 


-06] 


- 2.77866[- 


-06] 


8 


- 2.76786[ 


-07] 


+ 1.73880[- 


-07] 


9 


+ 1.24699[ 


-08] 


- 5.87474[- 


-09] 
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